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Many hydrated AB03-type perovskites are found to be 
fast proton conductors in the intermediate temperature range 
~200-500°C. Loading protons in the perovskite structure re- 
lies on acceptor-doping to the B site, whereby oxygen vacan- 
cies are created, which can be filled with hydroxyl groups in a 
humid atmosphere at elevated temperatures. 1 On a local scale 
the proton conduction process is composed of two elementary 
steps: (i) hydrogen-bond mediated proton transfer between 
adjacent oxygens, and (ii) rotational motion of the hydroxyl 
group in between such transfers. 1 Lhe long-range motion of 
the protons is a scries of such transfers and rotations, with an 
overall rate which depends on the local energy barriers. 

The long-range proton diffusion constant is commonly 
extracted from conductivity experiments via the Nernst- 
Einstein relation, but there is also one example of the use 
of pulsed-field gradient nuclear magnetic resonance (PFG- 
NMR) 2 . On microscopic length-scales, proton dynamics in 
hydrated perovskites have been investigated with frequency- 
resolved quasielastic neutron scattering (QENS), covering the 
ps time-scale, extended up to ~1 ns in some cases, in which 
mainly the local dynamical processes have been observed, 
although data have also been interpreted in terms of non- 
localized motions. 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9,10 Investigations of the proton dy- 
namics would, however, benefit greatly if one could extend 
the time-range to longer time-scales so that the microscopic 
diffusional process can be studied and information about the 
influence of the local structure and energy barriers on the pro- 
ton diffusion can be obtained. This is indeed possible by the 
use of neutron spin-echo (NSE), 11 ' 12 ' 13 with which the pro- 
ton dynamics can be studied over a microscopic length-scale 
in a wide time-range of ps-^is. However, to our knowledge, 
this technique has previously neither been applied to study 
the proton dynamics in hydrated perovskites nor in any other 
proton-conducting ceramic. 

In this Communication we demonstrate the applicability 
and potential of NSE to study proton dynamics in proton- 
conducting ceramics. This is exemplified by experiments per- 
formed on hydrated BaZro.90Yo.10O2.95 (10Y:BZO), a cubic 
perovskite with a relatively high proton conductivity. 1,14 The 
high proton conductivity together with a high thermodynamic 
stability make this material a promising candidate for use as 
electrolyte in intermediate temperature fuel cells. In addition 
to our experimental work, we use kinetic modeling based on 
first-principles calculations to assist the interpretation of the 
experimental data. 

The NSE experiment was performed at the IN 15 spectrom- 
eter at Institut Laue-Langevin (ILL) in Grenoble, France, 



with which a wide time-range over nearly three decades, ~0.2 
to 50 ns, was covered. The measured quantity in NSE is the 
polarization of the scattered neutrons as a function of the 
momentum transfer Q and Fourier time t, which is directly 
related to the intermediate scattering function I(Q,t), 1[i re- 
vealing the proton dynamics. For instance, for free diffusion 
the intermediate scattering function would show an exponen- 
tial relaxation, i.e. I(Q,t) = e~*/ T ^- ) , with the relaxation 
time r proportional to Q 2 , namely t -1 (Q) = DQ 2 where D 
is the diffusion coefficient, while a more complex motion would 
consequently result in a more complex shape of I(Q,t), see 
further below. 15 Data were obtained for the following temper- 
atures and Q-values: 473 K (Q = 0.3 A" 1 ), 521 K (Q = 0.3 
A- 1 ), 563 K (Q = 0.15 and 0.3 A" 1 ), 623 K (Q = 0.2 A^ 1 ) 
and 650 K (Q = 0.2 A" 1 ). 

Figure 1 shows the normalized I(Q,t), obtained in the ex- 
periment at 521-650 K. As can be seen directly in Figure 1, 
I(Q, t) is characterized by a single decay process in the time- 
range 1-5 ns for all four temperatures. On this time-scale we 
do not expect any other dynamics than proton motions, since 
the probed momentum transfers, Q — 0.15-0.3 A -1 , corre- 
sponds to a length-scale of ~20-40 A in real space. We can 
therefore safely assign the relaxation observed in the interme- 
diate scattering function to proton dynamics. Furthermore, 
we find that the decay in I(Q,t) is well described by a single 
exponential function for all temperatures, see Figure 1. For 
the temperature 563 K we have determined the relaxation rate 
for two Q-values and as can be seen in the inset in Figure 1 
the relaxation rate follows a Q 2 -behavior. The fact that the 
relaxation function in the experiment is well described by a 
single exponential decay, and is consistent with a Q 2 -behavior, 
shows that we are observing free translational diffusion of the 
protons. 15 This implies that already on a length-scale of 20 
A, which corresponds to a distance of ~5 unit cell lengths 
(the lattice parameter of 10Y:BZO is 4.20 A) 14 , potential lo- 
cal traps, or other "imperfections" in the structure that may 
severe the proton transport, have averaged out. Therefore, 
we have extracted a diffusion coefficient from the experimen- 
tally determined relaxation times assuming a r _1 (Q) = DQ 2 
dependence also for the other temperatures and the result 
from this analysis is shown in the Arrhenius plot in Figure 2. 
It is evident that the obtained values are consistent between 
the different temperatures, which provides support that our 
analysis is physically reasonable. 

To further corroborate our analysis of the NSE data we 
have also modeled I(Q,t) using a kinetic model based on 
first-principles calculations. 16 ' 17 The first-principles calcula- 
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FIG. 1: /(Q,t) at T = 650 (Q = 0.2 A^ 1 ), 623 (Q = 0.2 
A" 1 ), 563 (Q = 0.3 A" 1 ) and 521 K (Q = 0.3 A" 1 ). The 
solid lines are fits to I(Q,t) = e~*' T ' Q '. Inset in upper right 
panel: <2 2 -dependence of the relaxation rate at T = 563 K, 
where the line represents a fit to r _1 (Q) = DQ 2 with D = 
8.53- 10" 7 cmV 1 . 



tions were carried out within the framework of density func- 
tional theory (DFT); details on the settings can be found in 
ref. 16 and in the Supporting Information. From the DFT cal- 
culations the barriers to proton transfer and hydroxyl rotation 
in the vicinity of and far from dopants were assesed, and used 
as input parameters for the kinetic model. The calculated 
diffusion barriers far from Y-dopants are found to be 0.20 eV 
and 0.18 eV for the proton transfer and hydroxyl rotation 
motion, respectively. The binding energy to a Y-dopant is 
0.16 eV and we find that the influence of the Y-dopant on the 
energetics for the proton is quite extended in space, including 
both the first and second coordination shells. For the dop- 
ing level of 10% it implies that the region influenced by the 
dopants amounts to about half of the lattice sites and a con- 
ventional description in terms of a trapping diffusion model 
cannot be used. 

Figure 3 shows the calculated intermediate scattering func- 
tion, 7 ca ic. (Q, t), at T = 563 K and for a dopant concentration 
corresponding to 12.5%. Results are shown for momentum 
transfers Q = 0.3, 0.5, 2.0 A -1 and the long-range diffusion 
limit Q — > 0, in which the scattering function is given by a sin- 
gle exponential with a characteristic relaxation rate r _1 (Q) = 
DQ 2 , where D is the calculated diffusion coefficient. 15 Since 
the calculated scattering functions are plotted against Q 2 t 
they will collapse onto a single curve as long as we are in 
the long-range diffusion regime. As seen in Figure 3, this is 
clearly the case up to, at least, Q = 0.5 A^ 1 . Thus, for the 
Q- values probed in the NSE experiment we indeed expect to 
observe the long-range diffusional process. 

Only for larger Q- values J C aic.(Q,i) deviates clearly from 
the single exponential behavior. For example, at Q = 2.0 
A -1 in Figure 3 one can observe a two-step relaxation func- 
tion. This reflects a, so called, trapping behavior, i.e. the 
proton spends an extended time in the vicinity of dopant 
atoms before it diffuses further through the structure, thus de- 
creasing the long-range proton mobility. The view of dopant 
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FIG. 2: Diffusion constant obtained from NSE (•), first- 
principles calcualations (o), and conductivity measurements 
(line). The dashed line is an extrapolation of the latter result. 



atoms acting as well localized trapping centers has previ- 
ously been proposed based on quasielastic neutron scatter- 
ing data, 3 muon spin-relaxation experiments, 18 and computer 
simulations. 19 ' 20 ' 21 It should be noted that for lower dopant 
concentrations the transition from a single exponential to a 
more complicated relaxation will occur at smaller Q-values. 

From the kinetic model based on first-principles calcual- 
tions we have also determined the diffusion constant for five 
different temperatures. These values are inluded in Fig- 
ure 2 together with the values obtained from conductiv- 
ity measurements 22 for comparison. The conductivity mea- 
surements were performed at lower temperatures than our 
NSE measurements, but if one extrapolates the conductiv- 
ity data we find a good agreement between the proton diffu- 
sion constants obtained from the two different experimental 
techniques. The kinetic model on the other hand overesti- 
mates the diffusion constant by about one order of magnitude. 
However, the accuracy of the obtained potential energy sur- 
face is limited by the accuracy of the approximation for the 
exchange-correlation functional. It is known that generalized 
gradient approximations (used in the present study) have a 
tendency to underestimate energy barriers for proton transfer 
processes. 23 An underestimation of 0.1 eV, which is not at all 
unreasonable, would increase the rate by about a factor of 10 
at the present temperatures. Changes in the potential energy 
surface will also influence the value for the vibrational fre- 
quencies entering the expression for the prefactors for the in- 
dividual jumps. Furthermore, our modeling of the prefactors 
for the individual rates 16 can be improved by incorporating 
more vibrational degrees of freedom and the diffusion model 
we have employed assumes that consecutive proton jumps in 
the structure are uncorrelated. 

The fact that the diffusion constant obtained from the NSE 
experiment compares well to that derived from the conduc- 
tivity data is perhaps surprising since the two techniques 
probe considerably different length-scales. The diffusion con- 
stant probed in conductivity measurements reflects the bulk 
diffusion of protons between grain boundaries, with a length- 
scale of the order 1 /im, which is much longer than the length- 
scale over which the proton diffusion is investigated by NSE, 
~20-40 A in our case. This result suggests that already over 
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FIG. 3: Calculated intermediate scattering functions, 
^caic.(<2,£), for various momentum transfer values, Q, at 563 
K. 



a distance as short as five unit cell lengths we are measuring 
the true long-range proton self-diffusion in the material. Only 
for Q-values larger than 0.5 A -1 , corresponding to three unit 
cell lengths, trapping effects are expected in this material. By 
extending the Q-range to higher Q-values in the NSE exper- 
iment, it should be possible to observe the cross-over from 
a single exponential at low Q, typical for long-range proton 
diffusion, to a more complex form at larger Q-values, sug- 



gesting that several processes are reflected in I(Q,t). Such 
investigations may for example give information about the 
existence/nature of traps, an often debated subject for hy- 
drated perovskites. 

To conclude we present data from the first NSE experi- 
ment on a proton-conducting perovskite, namely hydrated 
BaZro.90Yo.10O2.95, and demonstrate that NSE spectroscopy 
can be beneficially applied to investigate the protonic 
self-diffusion in this class of materials, and most likely also 
in many other types of solid-state proton conductors, such 
as proton-conducting polymers, 24 solid acids, 25 and alkali 
thio-hydroxogermanates. 26 Furthermore, we show that NSE, 
which indeed is the only neutron scattering technique that 
gives access to the long time-scales needed to accurately 
investigate the comparatively slow proton diffusion in hy- 
drated perovskites, can be powerfully combined with kinetic 
modeling based on first-principles calculations, since both 
techniques cover the same time- and length-scales. 

Acknowledgment. We acknowledge financial support 
from the Swedish agencies VR, NFSM, and SSF via the 
ATOMICS program. Allocations of beam time at ILL, and 
of computer resources through the SNAC are also gratefully 
acknowledged. 

Supporting Information Available: Details of the sam- 
ple preparation, the NSE experiment, and the kinetic mod- 
eling based on first-principles calculations. This is available 
free of charge via Internet at http://arXiv.org/. 



1 K. D. Kreuer, Annu. Rev. Mater. Res. 33, 333 (2003). 

2 K. D. Kreuer, T. Dippel, Y. M. Baikov, and J. Maier, Solid 
State Ionics 86-88, 613 (1996). 

3 R. Hempelmann, C. Karmonik, T. Matzke, M. Cappado- 
nia, U. Stimming, T. Springer, and M. A. Adams, Solid 
State Ionics 77, 152 (1995). 

4 T. Matzke, U. Stimming, C. Karmonik, M. Soetramo, 
R. Hempelmann, and F. Guthoff, Solid State Ionics 86- 
88, 621 (1996). 

5 B. Grofi, C. Beck, F. Meyer, T. Krajewski, R. Hempel- 
mann, and H. Altgeld, Solid State Ionics 145, 325 (2001). 

6 C. Karmonik, R. Hempelmann, J. Cook, and F. Guthoff, 
Ionics 2, 69 (1996). 

7 M. Pionke, T. Mono, W. Schweika, T. Springer, and 
H. Schober, Solid State Ionics 97, 497 (1997). 

8 D. Wilmer, T. Seydel, and K. D. Kreuer, Mater. Res. Soc. 
Proc. 972, 15 (2007). 

9 M. Karlsson, A. Matic, D. Engberg, M. E. Bjorketun, 
M. M. Koza, I. Ahmed, G. Wahnstrom, P. Berastegui, 
L. Borjesson, and S. G. Eriksson, Solid State Ion. 180, 
22 (2009). 

10 A. Braun, S. Duval, P. Ried, J. Embs, F. Juranyi, 
T. Strassle, U. Stimming, R. Hempelmann, P. Holtappels, 
and T. Graule, J. Appl. Electrochem. doi:10. 1007/210800- 
008-9667-3 (2008). 

11 F. Mezei, Neutron Spin Echo: Lecture Notes in Physics, 
vol. 28 (Springer, Heidelberg, 1980). 

12 F. Mezei, Zeitschrift fur Physik 255, 146 (1972). 

13 F. Mezei, C. Pappas, and T. Gutberlet, Neutron Spin Echo 



Spectroscopy: Basics, Trends and Applications (Springer, 
2003). 

T. Schober and H. G. Bohn, Solid State Ionics 127, 351 
(2000). 

M. Bee, Quasielastic Neutron Scattering (IOP Publishing, 
Bristol, D.J. Millen ed., 1988). 

M. E. Bjorketun, P. G. Sundell, and G. Wahnstrom, Fara- 
day Discuss. 134, 247 (2007). 

M. E. Bjorketun, P. G. Sundell, G. Wahnstrom, and D. En- 
gberg, Solid State Ionics 176, 3035 (2005). 
R. Hempelmann, M. Soetratmo, O. Hartmann, and 
R. Wappling, Solid State Ionics 107, 269 (1998). 
R. A. Davies, M. S. Islam, and J. D. Gale, Solid State 
Ionics 126, 323 (1999). 

M. S. Islam, P. R. Slater, J. R. Tolchard, and T. Dinges, 
Dalton Trans, pp. 3061-3066 (2004). 

M. S. Islam, R. A. Davies, and J. D. Gale, Chem. Mater. 
13, 2049 (2001). 

K. D. Kreuer, S. Adams, W. Munch, A. Fuchs, U. Klock, 

and J. Maier, Solid State Ionics 145, 295 (2001). 

V. Barone and C. Adamo, J. Chem. Phys. 105, 11007 

(1996). 

M. F. H. Schuster and W. H. Meyer, Annu. Rev. Mater. 
Res. 33, 233 (2003). 

S. M. Haile, D. A. Boysen, C. R. I. Chisholm, and R. B. 
Merle, Nature 410, 910 (2001). 

M. Karlsson, A. Matic, I. Panas, D. T. Bowron, S. W. 
Martin, C. R. Nelson, C. A. Martindale, A. Hall, and 
L. Borjesson, Chem. Mater. 20, 6014 (2008). 



Using neutron spin-echo to investigate proton dynamics in proton-conducting 

perovskites - Supporting Information 



Maths Karlsson, 1 ' 2 Dennis Engberg, 2 Marten E. Bjorketun, 2 Aleksandar Matic, 2 Goran Wahnstrom, 2 Per G. 
Sundcll, 2 Pedro Berastcgui, 3 Istaq Ahmed, 4 Peter Falus, 5 Bela Farago, 5 Lars Borjesson, 2 and Sten Eriksson 4 

European Spallation Source Scandinavia, Lund University, SE-221 00 Lund, Sweden 
2 Department of Applied Physics, Chalmers University of Technology, SE-412 96 Goteborg, Sweden 
'Department of Inorganic Chemistry, Arrhenius Laboratory, 
Stockholm University, SE-106 91 Stockholm, Sweden 
4 Department of Chemical and Biological Engineering, 
Chalmers University of Technology, SE-412 96 Goteborg, Sweden 
5 Institut Laue-Langevin, 6 rue Jules Horowitz, BP 156, 38042, Grenoble Cedex 9, France 

(Dated: June 12, 2009) 

Sample preparation: The BaZro.90Yo.10O2.95 sample was prepared by mixing stoichiometric amounts of BaC03, 
Zr02 and Y2O3. The oxides were heated to 250°C overnight to remove moisture prior to weighing. Milling was 
performed manually using an agate mortar and a pestle. The finely ground mixtures were fired at 1000°C for 8 
hours and subsequently ground and pelletized using a 13 mm diameter die under a pressure of 8 tons. The pellets 
were sintered at 1200°C in air for 72 hours. After sintering, the pellets were reground, compacted, and retired at 
1500°C for 48 h. Finally, the pellets were finely reground to powders. The charging with protons was performed by 
annealing the powder samples at 300°C under a flow (12 ml/min) of Ar saturated with water vapor at 76°C for 
10 days. X-ray diffraction measurements, performed at ambient temperature on a Siemens D5000 powder diffrac- 
tometer (CuKai = 1.5406 A), revealed a cubic structure of the hydrated material, in accordance with the literature [1]. 

NSE experiment: The NSE experiment was performed at the IN15 spectrometer at Institut Laue-Langevin 
(ILL) in Grenoble, France. We used neutrons with wavelengths centered at 10 A. The instrumental resolutions were 
calibrated using a graphoil elastic scatterer. The data analysis was performed using the standard IN15 data package 
and FFJDA [2]. The powder sample was loaded in a vacuum tight Al container, which was coated with a 100 nm thick 
layer of Pt on the inside to avoid corrosion. The sample thickness was chosen to 6 mm, to obtain a total scattering 
of neutrons of approximately 10%. It should be noted that the experimental data is customarily normalized through 
a division with I(Q,t = 0) [3], but in our case this value was close to zero, so normalized scattering functions have 
instead been obtained from fits of the unnormalized data, I cxp .(Q, t), according to I eKp .(Q, t) — a(Q) + b(Q) ■ I(Q, t), 
where a and b are Q-dependent fitting parameters. For each spectrum, the data was logarithmically summed into 8 
data points, in order to increase the statistics. 

Computational details: The first-principles calculations were carried out within the framework of density func- 
tional theory, based on the plane-wave pseudopotcntial approach as implemented in the Vienna ab-initio simulation 
package (VASP) [4,5]. The electron-ion interactions were described by the projector augmented wave method [6]. For 
the exchange-correlation part, we used a generalized gradient approximation (GGA) due to Wang and Perdew [7]. 
All calculations were performed non-spin-polarizcd with a plane- wave cutoff of 400 eV. 

In most energy barrier calculations, the hydrated Y-doped BaZr03 sample was represented by a supercell consisting 
of 3x3x3 primitive BaZr0 3 units plus one Y ion, which corresponds to a doping level of 3.7%, and the Brillouin 
zone sampling was performed using a 2x2x2 k-point grid. In order to ensure the transferability of the results to 
higher dopant concentrations, similar to that of the "experimental" sample, a few key barriers were also calculated 
using a supercell with 12.5% doping. Since only minor changes in the barrier heights were observed upon increasing 
the doping level, we feel confident that reliable predictions about 10Y:BZO can be made based on the 3.7% data. 
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